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Abstract 

The investigation of nonlinear dynamical systems of the type 

X = P{x, y, z), y = Q{x, y, z), z = R{x, y, z) 

by means of reduction to some ordinary differential equations of the second 
order in the form 

y" + ai(x, y)y'^ + 3a2(x, + 3a3(x, y)y' + ai{x, y) = 

is done. The main backbone of this investigation was provided by the theory 
of invariants developed by S. Lie, R. Liouville and A. Tresse at the end of the 
19*^ century and the projective geometry of E. Cartan. In our work two, in 
some sense supplementary, systems are considered. The first one is the Lorenz 
system 

X = a{y — x),y = rx — y — zx , z = xy — bz 
where a, r, b are parameters and the second one is the RoBler system 

X = —y — z,y = X + ay , z = b + xz — cz 

where a, b and c are parameters. The invarinats for the ordinary differential 
equations, which correspond to the systems mentioned abouve, are evaluated. 
The connection of values of the invariants with characteristics of dynamical 
systems is established. 
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1 Introduction 



The first broad investigation of invariants of the second order differential equation 

y" = f{x, y, y') (l) 

under some general point transformations 

y = V{u,v), (2) 

where and rj are arbitrary smooth functions was done in the work of R. Liouville . 
He found some series of absolute and semi-invariants and discovered one procedure 
to build other invariants of higher weights. A further consideration of the same 
problem from the point of view of infinitesimal transformations was introduced by 
S. Lie 0, 1^ and completed by his student A. Tresse for second order differential 
equations (|I|) with arbitrary smooth functions f{x,y,y') in 0]. The geometrical 
description of these results was given by E. Cartan when he introduced the 

new idea of projective connections. It was important for the theory of invariants 
to investigate the invariants themselves and their possible connections among each 
other, because it enables us to say whether our equation (|1]) allows some infinitesimal 
group of point transformations or not @]. In projective geometry it is important to 
know if our invariant is equal to zero or not, because this circumstance is responsible 
for the existence of corresponding geometrical properties. 

We will use these ideas to investigate the nonlinear dynamical system of the type 

X = P{x,y,z) , 
y = Qix,y,z) , 

z = R{x,y,z). (3) 

At first we conduct some "projections" of this system on, for instance, the {x, y)- 
plane or the {y, 2;)-plane or some nonplanar surface. After that we deal with some 
ordinary differential equations. 

In all our investigated cases these equations have the form 

y" + ai{x, y)y'^ + 3a2(x, y)y''^ + 3a3(x, y)y' + 04 (x, = . (4) 

It is well known that this form will be keeped under the point transformations 
and invariants of such equations can be easy investigated. Hereforth we can look for 
conditions for the parameters in in which one or more invariants are equal to zero. 
We hope that exactly these conditions will be able to show us where integrable cases 
lie. It is well known that some typical Riemann surfaces cannot yield a satisfactory 
corresponding geometrical interpretation of the equation (^) in general, because 
metrical properties are noninvariant under the point transformations. The adequate 
geometrical image is only possible with the help of the notion of normal projective 
connection. 
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According to E. Cartan ^ we will use the following notations. For each second 
order differential equation we will look at the ccorresponding three dimensional 
vareity M^{x,y,y') of elements of first order with local coordinates {x,y,y') and 
projective connections. The structure of this variety can be described through the 
1-form components of projective connection 

w\ w\ - wl, wl i ^ J, i,j = 1, 2, 3 

and 2-forms of curvature and torsion 

However, in the case of equations of the type (^) a torsion free variety with normal 
projective connection is produced and our variety M^{x, y, y') can be represented as 
a direct product 

M^{x,y,y') = V\x,y)xSK 
The 2-forms of curvature and torsion satisfy the following conditions 



= Q^ = 0, Ql = Ql = 0, Q\-Q^o = Ql-Q°^ = 0, 



(5) 



where 



o^^«2 S^as d^di da^ dai da^ dai 

= 2—— - + 2ai^ + 3ai— Sag— — 

oxoy oy^ ox^ oy oy ox ox 

das da2 , / d'^a^ d'^a2 d^a^ dai 
oa2^ h ba2^ h 2-— —5 — ^ — 204 



h = 



dy dx \ dxdy dx"^ dy"^ dx 

da^ da2 da^ da2 da3\ . 

ai— + 3ai— + 3a2— + 3a3— 603— y, (6) 

dx dy dy dx dy J 

db 
dy' ■ 



The components of the normal projective connections also have a very simple 
form 

= dx, w"^ = dy, wf = a^dx + a^dy, wl = —a2dx — aidy, 
wl = —wl = —a^dx — a2dy, w^ = Il^idx + n?2d|/, W2 = Il^ida; + Il22dy, 

with 

^11 = 2(03 — a2a^) + a^x — o-Ay ? 
■^22 ~ 2(a2 — difls) + aix — a2y , 
ni2 = = 0203 ~ aia^^ + a2x — cisy ■ 
The Cartan's structure equations have the form 

dw* = w^ Aw\, 

dwj = w^ A w^ + w^ A wl — djw'' A wl , 

dw° = w'' Awl - I^m'^'' j, k,l = 1,2. 



The tensor of projective curvature has two components i?5'i2 -R2125 ^^Sit will be 
noted as —Li and —L2 respectively: 



dy dx 



L2 = + ^ - «in?i - a,Ul, + 2a,U% . (7) 

Both these values were discovered by R. Liouville much earlier than by E. Cartan 
and used for the construction of the semi invariant z/5 with the weight equal to 5. 
This invariant is a semi invariant. This means that after some point transformation 
(g) the new 1/5 is equal to 

z75 = AV5 , 

where A 7^ is the functional determinant of our transformation in opposite to an 
absolute invariant where we have f = r. The degree of this determinant is exactly 
the weight of the corresponding invariant. The value of can be calculated by 
means of the following formula: 

= L2{LiL2x — L2Li^)-\-Li[L2Liy — LiL2y)—aiL\ + ?)a2L\L2 — 'ia^LiL\-\-a4^L\. (8) 

If 1^5 7^ we can use the recursive formula of R. Liouville [|l| to build the series 
of semi invariants with higher wights: 

Urn+2 = Li-^ - L2-^ + mUm{L2x - Liy), 171 > 5 , (9) 

and after that construct the absolute invariants 

tm = l^m^^"^'^ ■ (10) 

In the case z/5 = 0, R. Liouville has found other seires of semi invariants. At first, 
he found some invariants of weight 1: 

X 

wi = Yl [-^2("'-^2 - OiLi) - R2{Ll)y + L2R2y - L2R2{aiLi - 02-^2)] 
for L2 ^ 0, and 

w^i = 4 [Llia'L, - a"L2) + - + URi{a^L^ - 04^2)] (11) 

for Li 7^ 0, where 

Ri = L1L2X ~ L2L1X + (I2LI — 2a^LiL2 + (34-^2, 

R2 — LiL2y — L2Liy + (1iL\ — 202-^1/72 + Q,3-Z^25 

a = a2y — aix + 2(0103 — al), a' = a^y — a2x + «i«4 — 0203, 
a" = a^y — a^x + 2(0204 — O3) . (12) 
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Also in this case there exists one recursive formula to construct the series of semi 
invariants with higher weights 

Wm+2 = Li-^ - L2-^ + mWm{L2x - Liy) (13) 

and the corresponding absolute invariants read 

Wm+2 ^ . ,x 

Um+2 = ^71+2 ' "2 > 1 . (14) 

The case that both projective curvatures Li and L2 are equal to zero is non inter- 
esting, while the equation (|) can then be reduced to y" = . If also Wi = there 
are other series of invariants. The first one has weight equal to 2 and is given by 

. _3Ri^dL2 _dLa 
Li dx dy 

and other semi invariants of this series by 

di2m J di2m , ^ . /^Lg dLi\ 

%2m+2 = Li— L2— h 2m«2m TT- ] , [m > 1). (16) 



oy ox \ ox oy 

The corresponding absolute invariants can be found through the relation 

j2m = «2m^2 • (17) 

If for our equation (^ the semi invariant 

z/5 = (18) 

and Li, L2 that we can immediately construct the first integral of this 
equation 

y' = -^. (19) 

For many physical systems like (B|) those partial solutions are very interesting. 



On the other hand side the condition (pTSD has also a deep geometrical sense. 



Namely, if and only if the condition (|T8D is fulfilled it will be possible to immerse 
our variety V'^{x,y) of the normal projective connection into the real projective 
space P^(M). In all other cases the immersion will be possible only in P^(]R) [^]. 

The immersion into P^(]R) will be realized on one of the simplest surfaces - on 
developed surfaces. It is well known that the coresponding surfaces for the equation 



y" = is the projective plane. In this sense the equation (^ with the condition ([18|) 
are really the next simple ones. As we mentioned above, the condition that some 
invariant is equal to zero must have some geometrical sense, is in the cases Wi = 
and ^2 = not known to us. 

The first investigation of dynamical systems by the described method was done 
in 11,11]. In the present work we deal with different aspects of the Lorenz system as 
well as with the RoBler system. 
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2 The Lorenz system 



In the work ||T5| Lorenz discussed the problem of the representation of a forced 
dissipative hydrodynamic flow and suggested the following model for the description 
of that system: 



dx 

dt 
dy 

dt 
dz 

dt 



a{y - x) , 
rx — y — zx , 

xy-bz, (20) 



where a, r, b are some parameters characterizing the flow. Now, this system is well 
investigated by some analytical and numerical methods. It posseses several dynam- 
ical states corresponding to different regions in the parameter space. However, up 
to now no one analytical criterion is known which allows to decide wheter the for a 
given set of parameters the solution will be a regular or a stochastic one. 

The phase space of (^) is three dimensional and it is separeted due to the 
integral curves of the following two equations 

dy rx — y — xz dz xy — bz 
dx a{y — x) ' dx a{y — z) 

After eliminating one of the functions, z, for instance, and introducing new notations 
we get an equivalent second order differential equation of the type (^ 

d^y Sdy'^ f 1\ 4 , V 3 n 3 4 n 2 3 n lnn\ 
+ \ay — h exy H 7?/ — px y — px y = (22) 



dx^ y dx \ X J dx x 
with 

b + a + 1 1 b{a + l) a + 1 b{r - 1) 

« = , P = 7 = 2 — ' " = ' ^ = 2 — • (^23) 



Roughly speaking, this is the projection of the system (^OD on the (x, ?/)-plane. The 
equation (^) was investigated in the paper pj. The semi invariant z/5 has the form 

1^5 = UiZ + U2't'^ + W4 (24) 

with z = x'^, t = l/{xy) and 

Ui = aP{10a — c? — 65), 
2 

^2 = -a(a + 3(5)(2a — 35), 
9 

= «( — (2^2-97) -2e(2a- 3(5)). (25) 
9 
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The projective curvatures are 

2 1 

Li = 3(ex - (3x^)y^ - -a% + —(2a - 35) , 

3 3a; 

L2 = -. (26) 

y 

The invariant 1/5 is equal to zero in two cases. First, for a = 0, i.e., a = — 6 — 1. 
This case is nonint cresting because it cannot be reached for physical values of the 
parameters. In the second case we have = for arbitrary r and a = 2h — 1 with 
6 = and h = 2/3. This is a well known case of regular behaviour of the system 

In the following we shall continue the investigation of and calculate the other 
invariants as well as the invariants of other projections of the system (pO|). 

R. Liouville [|l| found that in the case z/5 = and Li, L2 7^ we can construct 
the first integral of the equation (^): 

y' = -^. (2T) 

Consequently, for the projection on the (x, ?/)-plane with the parameters 

a = -1/5, b = -16/5, r = -7/5 (28) 

we get the equation 

y' = -^(4 - 32xy + 192x^y^ - 25a;V) , (29) 

X 

where y{x) is the solution of equation (|22|) . Let xy = s, then we get Abel's equation 
of the first kind [16] 

s. = f2ix)s^ + Ux)s^ (30) 

and after the substitution 

s{x) = w{x)r]{0 - M"^) = ^' M^) = ~^~ X^' 

(X^ \ 63 

^ = y /3(x)w;^(x)dx = (l68 + 25x^)^xi - y y"x^(168 + 25x^)-idx 
we have the following equation for the function ?7(^) 

with 

, 256 256 + 225x2 
I(x) = 



27 (168 + 25x2)47/21 
7 



and ^ as in equation (|3lD. 

We believe that this is a new first order integral of the Lorenz system, although 
it may be of purely mathematical interest due to the choice of parameters given by 
equation (^). 

Another case when a first integral of type ( P7D exists corresponds to the regular 
behaviour of the Lorenz system, that is to 

a = 1/3, b = 2/3, r arbitrary. 

Then we have 

9 

y' = 1/^(4 + 3(1 - r)xy + -x^y) . 



This is the well known Abel equation (|30D with 

/3 = 3(l-r)a; + ^x^ = 4 

and after a substitution in analogy with (^ip we get the canonical expression 

4 

y{x) = w{x)r]{^) - 



3(3(1 -r)x + 
-16 

3(3(1 fx 



^ = exp ( / X . Q 3x dx 



with 



and 



V(0 = + H^) 



To our knowledge, there are no such solutions in the literature. 

Let z/5 7^ 0, then our two dimensional variety V2{x,y) can be immersed into 
P^(M) only. On the other hand side we can calculate the invariants of higher weights 
and look for conditions on the parameters that one of that invariants vanish or 
all invariants are functions of one of them. In the first case we can hope that our 
immersed surface has some special geometrical properties. In the second case we can 
decide wheter our differential equation (^) possesses some group of infinitesimal 



transformations or not. 

We calculated 1/7, z/g, z/n, and some resultants of z/5 and 1/7 — -R57, z/5 and uq — R^q, 
uj and z/g — -R7g, as well as the resultants of i?57, -Rsg, i?7g. The results of these 
calculations are quite voluminous and only some first of them will be listed here. 
This gives us enought insight into the structure of them. 

The semi invariant of weight 7 for (|2^) equals to 

Z/7 = (f 1-2^ + V2Zt^ + V^Zt + f 42; + t>5t^ + V^t^ + Vjt^ + t>gt + t>io) / 1 (32) 
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with z = x^, t = l/{xy) and 

vi = 30/3-ui, t>2 = 2a/5(a; — 2)(a^ + 185), = —a'^ui, 

Vi = 10ap{ — {2a^-9-f)+3e{a'^-Ua + 126)),vr, = -{a + 36)u2, (33) 
3 9 

We = — + 36){2a — 36), vj = — — ea'^fg, vg = —a^Ui, fio = 30e'U4 

with ui,U2,U4 defined in (^Sf) . For the semi invariant of weight 9 we found the 
representation 

Ug = {wiZ^ + W2Z^t^ + Wj.z'^t + Wis^z'^ + W^zt'^ + W^zt^ + WTZt^ + Wg-Zt + 

WgZ + Wiot® + Wllt^ + Wi2t'^ + Wigt^ + Wut^ + WiqI + ^17)/^^ , (34) 

with z and t as above and 

wi = mop^i, W2 = l0aP^{-U8a'^ + 22a^ + 210a + S + 3a^5 - 3065^), 
W3 = 260a^/5Mi, W4 = 130a/3^(a2(2a2-97) + 18e(a2-12a + 95)), 

^5 = '^a'^(3{98-15a) + 2a/3{{a^ delta{2 + 6) +6^{78a- 1386)), 
9 

We = 52a^/?(a2(2 + a) + 6(5(lla - 245))/9, 

?2;7 = 10a/3(-a^(— - 14a + 105) +4a2e(ii _4a) + 
9 

35e(20452 - 104a - a^) + a%(7a + 35)), 

= 260a'^l3{-a\2a^ - 97) + 9e{a^ - Ua + 125)), (35) 
9 

wg = -130af3e{2a^{2a'^ - 9j) + 9e{a'^ - 18a + 185)), 

2 68 
Wio = -(a + 35)^-U2, Wii = — a^(a + 35)'U2, W12 = 2(14a'^ — 23e(a + 35))'U2, 
9 27 

Wi3 = — 416ea^'U2, 

14 

wi4 = 10a(— a^(2a2-97)-2ea^(3a- 55) + 34e2(2a2 + 3a5-952)- 
81 

a^e7(a + 35)), 
WiQ = —2Q0a^eu4, Wn = 1170e'^M4 



with Ui,U2,U4 from (pSl). 

The first two of the resultants -R57 and R^g have a short form. Using the notions 
from (^) and (|33| - p5|) we have 

-R57 = l^5ulvi + ulvi + ulvio - U1U4V4 + Z/5(2u3£t4fi - U1U3V4) + t'^{u\vi - U1U2V2 + 
u\v^) + t^{u\v(i - U1U2V3) + t^{2U2U4Vl - U1U4V2) + V^{2U2U3VI - 
U1U3V2) - MlM2'i^4 + U1V7 + t{iyjulvs + ulvg - V^UiU^V^ - U1U4V3), 
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ulwio - Uiulw2 + u\u2W^) + t^{-u\wii - Uiu\w':i + ulu2We) + 

t yiu^U^Wi — U\W\2 — 2U1U2U4W2 — + + Z>'5(3-U2W3Wl — 

+ ulu'iW^) + u\u2W'j) + t^(-nitf;i3 - 2-Ui'U2^t4^l'3 + S^£t4W6 + 
V^{-2uiU2U'iW'i + ulu^We) + u\u2Ws) + t(-M!^Wi6 - Z/5M1M3W3 - 
■U1M4W3 + ■U^M4W8 + i^5(-2'Ul'U3M4'^'3 + ■U^tt3W8)) + uluiWg + 
t {3U2U^Wi — U^WiA — UgUiWib — + Z/5(3U2M3Wl — — 

2'Ui'U2M4W4 + ulu^Wj + Z/5(6-U2'S3'U4Wl — 2uiU^U4W2 — 2U1U2U3W4, + 
ulusWj) + ulu2Wg) + ^^{SUsuIwi — 2UiU2,U4W4 + ulu^Wg) . 



Based on the analysis of the formulas ( p2ip5D we raise the supposition that each 
of the following higher invariants z/7, z/g, Uu can be equal to zero for the same values 
of the parameters as z/5. We failed to find different conditions on the parameters of 



^5 = 
z/7 = 
z/g = 



the system (pOl). Because all semi invariants of this series can be found by use of the 
recusive relation (|^) this supposition will we certainly true for all of them. Exactly 
the same can be said about the structure of the resultants. Let us remark that all 
these formulas become considerably shorter and simpler in the case a = 2b — 1 and 
only in this case we obtained a transparent relation between the invariants. Really 
we see, that for a = 26 — 1 

27b\-2 + 3b)z 

(26-1)5 ' 
16262(-2 + 3b)z{5b - 5bq + 5bH + bt^ - 2bH'^ + 5z) 

(26-1)7 ' 

106V(3(1 - 26)(1 -q) + 76^) + 2663^^(1 - 26)+ 
bh\l - 26)^ + 1306z(l - g) + 1306^^2 + 406^^2(1 - 26) + 65z'^)) , 

^^|^p^(10z/,2(-l + 26)^ + 54z/56^(2 - 36)(-5(l - q) - 

5bt - (1 - 2b)t^) + 9z/7(2 - 36)6^(1 - 2b)h)) . (36) 

Also, in this case we have a good correspondence with the physical interpretation 
of the system. The case a = 26 — 1 is the well known integrabel case and the system 
does not possess any stochastic behaviour. 

The investigated projection at the (x, ?/)-plane was the simplest one. In the fol- 
lowing we shall take more complicated projections and look for conditions on the 
parameters. In the first out from the more complicated cases we start with the sub- 
stitution ^ = — 62 in the system (|2l|) and eliminate the variable y from both 
equations. In this way we obtain the (x,^)- projection. For convenience we shall 



R57 — 
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use the notation y instead of ^. The second order differential equation which we 
obtain in this way has the canonical form given by equation with functions 
ai(x, y), i = 1, 4 as follows: 

cr(r — l)x ^ ax ax^ 
' h^y"^ b^y b^y"^ ' 

^ _ a+1 {2a-b){r- l)x^ {2a - b)x^ ^ {2a - b)x^ 1 

' 3by b'^y'^ b^y b^y'^ 3y ' 

, , {r -3){2a -bfx^ 2{2a - b){a + l)x 

(^3{x,y) = 7^ ^ rr + 

b^ay^ Sbay 
{2a - 6) V {2a - bfx^ ^ {2a - b)x 
b^ay b^ay"^ 3ay ' 

a {x y) - ~ ^)^^^ _ (2^ - b)^x^ ^ (a + l)(2a - 6) V _ 
' b^a'^y'^ b^a'^y a'^by 

(r - 1)(2(T - 6)V 2a -b 

In this projection the projective curvatures have the form 

Li = ^,"1 ~ i^x%{b - 2af{l - r) + Qx\b - 2af + 2x^yb{-b{l + bf+ 
3b^a''y^ 

2a{l + 86 - 6r{b - a) - Aa + a^ + 3ba)) + 30x^ya{b - 2a) - 
y%'^a{l + b + a)- 3x'^y^a{b - 8a)) , 

L2 = —^\-9x%{b-2a)\l-r)-9x\b-2a)^ + yb{2b + b^-b^+ 
3b'^ay'* 

a{-A - 106 + 36^ + 16a + Aa^ + 12r{b - 2a))) - 30x'^ya{b - 2a) + 
3y^a{b - 8a)] , 
and the semi invariant z/5 is given by 
(l + 6 + a)(2a-6) 
276V52/10 

(2 + 6 - 6^ _ + 6(T + 9r(T + 2a'^)x^ + 9b{2a - bf{3b + + 36^ + 
a - 346(7 - 176V + 186ra + 20(7^ - 196(T^ - ISra^ + + 
8\a{a - 2b){2a - 6) V + 186a(2a - 6)^(-36 - 66^ - 36^ - a + 346(t + 
176V - \8hra - 20a^ + 19ba^ + 18ra^ - a^)x^y + 
324(26 - a)a'^{2a - b)^x'^y + b^a^{2a - b)^{-2 -b + b^ + ba -ba - 
9ra - 2(7^) (3 + 26 - 6^ - 10a + 2ba + 16ra + 3a'^)xy'^ + ba'^{2a - 6)^ 
(336 + 486^ + 156^ + b9a - 2066(t - 1396V + 1446ra - 170^^ - 9bba'^ + 
288ra2 + h9a^)x^y'^ + I8a'^{2a - 6)^(36^ - 446a + lOa^)^^?/^ + 
26a^(2a - 6) (76 + 26^ - 56^ - 59a - 446a + 396V + 366ra + 170a2 
- 156a2 - 288ra2 - 59a^)xy^ + 36(6 - 8a) (6 - 2a)a^(6 + 2a)xV - 

9a^(8a - h)^xy' 



^5 



362(2a - 6)^(-l + 6 + 26^ - 5a + 6a + 3ra - a^ 
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In this case we can succeed in beeing identically equal to zero for 



a = -1-b and b = 2a (38) 

only. The first condition belongs to the mentioned above unphysical case, the sec- 
ond one represents the integrable case. The surface V2{x,y) can be immersed as a 
developeable surface in the projective space P^(R). 

For b = 2a we can investigate other series of invariants. The following ones 
(see above (p!lD -(|T3|)) with the first weight equal to 1 will be absent too because of 
Wi = for the considered values of the parameters. The next series of semi invariants 
i2k, k = 1,2, ... (p!5|)-([T6|) reads 

. _ 3 • _ 9 • _ 27 . 243 



4crV' 8(tV' 8(tV' 16(tV' 

The projective curvatures have a very simple form in this case 

3x 



Li = 0, L 



2 



Aa'^y'^ 

and the functions Oi (|37|) are given by 



oc 1 — (y 

Qi = o 9 o - + 2gcr - 2a), a2 = — , 03 = 04 = 0, 

8cr^?/"^ bay 



The corresponding absolute invariants j2k (ITtD are constants: 

j4 = 2, is = 8, J8 = 48, .... 

The last of the investigated projections of the Lorenz system is the (y,z)- projec- 
tion. The coefficients of the ordinary differential equation (^) have in this case the 
following form: 



ai{x,y) 



a2{x,y) 



,2^3 



ay — by X — bay x + b rx — b x 

(— 1/2 -|- brx — bx'^Y 
y{by'^ + ay"^ - Sray'^ - b'^rx + 2brax + 3ay'^x + - 2bax^) 



3(— ?/2 + brx — bx'^Y ' 
a^{x, y) = (ry^(l — b — 2a + 3ra) + 6r^x(— 1 + b — a) + y'^x{l + 2a — 6ra) + 

rx'^{2b - 6^ + 2ba) + 3ay^x^ - bx^{l + a))/{3{-y'^ + brx - bx^f) 
y{r'^a — r^a + y'^ — brx — 2rax + 3r'^ax + bx"^ + ax"^ — 3rax'^ + ax^) 



a4{x,y) 



_y2 _|_ Jjir-j, _ bx'^y 



For convenience we substituted z for x. This is one of the complicatest cases and 
all formulas are quite large. Therefore we list only the projective curvatures Li and 
L2 her. The first projective curvature has the form: 

Li = ?/(rV(4 -36^ + 6^ + 20a- 266ct + 86V- 12r(T + 246rfT- 96V(t + 4(t2- 



12 



2ba^) +y\2- 26^ + 2a - 2ba - 3ra + 6bra) + r^x{-Ab + 3b^ - 6^ + 
16ba - 106V + 6V - 2Abra + Ub'^ra - Aba^ + 26^^) - ry^x{8 + 46 + 
106^ - 26^ - 40a + 706a - 296V + 36ra - 966ra + 396Va - Sa^ + 66a^) + 
3ay^x{l - 26) + rV2(126 - 66^ - 136^ + 56^ - 486a + 466V - 66V + 
966ra - 576Va + 126a2 - 86V2) + iy^x'^^l - 6 - 6^ + 6^ + 5a - 116a + 
66V - 9ra + 306ra - 156Va + a^ - ba^) + brx^ {-12 + 126 + 146^ - 106^ + 
48a - 626a + 156V - 144ra + 996ra - 12a2 + 106a2 + ayV^(12 - 486 + 
306^) + x^(46 - 66^ - 46^ + 66^ - 166a + 266V - 106V + 966ra - 756Va + 
46a2 - 46V2) + 36ax^(-8 + 7b))/{3{-y^ + brx - bx^), 

the second one 



L2 = ry\6 - 6 - 56^ + 26^ - 10a + 156a - 56V + 12ra - 246ra + 96Va - Aa^ + 
2b(r^) + 3ay^(l - 26) + 2br^y^x{-2 + 6 + 26^ - 6^ - 12a + 106a - 26V + 
12ra - 66ra + 2(1^ - ba^) + 2y^x{-3 + 26 + 36^ - 26^ + 5a - 136a + 86V - 
12ra + 366ra - 156Va + 2(r^ - 2ba'^) + 6W^(-2 - 6 + 6^ - 2a + 6a) + 
bry^x'^{8 - Qb - 86^ + 66^ + 48a - 646a + 156V - 72ra + 456ra - Sa"^ + 
6a?/V(2 - 66 + 56^) + 6W^(6 + 26 - 46^ + 6a - 46a) + 2by^x^{-2 + 26 + 
26^ - 26^ - 12a + 226a - 106V + 36ra - 176ra + 2a^ - 2ba^) + 
6Vx^(-6 - 6 + 56^ - 6a + 56a) + 6ayV'^(-24 + 216) + 
26V^(1 - 6^ + a - 6a))/(3(-y2 + brx - bx^)- 

It is easy to observe that both curvatures are much simpler by a = and are both 
equal to zero in case 

a = 0, 6 = — 1, r arbitrary. (39) 

That means that the corresponding second order differential equation, the {y, z)- 
projection of the system (0), are equivalent to the simplest equation y" = and 
can be integrated. 

3 The Rofiler System 

For the Rofiler system [pi|-pi 



X = -y - z, 
y = x + ay, 
z = b + xz — cz , 

a, 6, c being parameters, we investigated 3 projections. In all of them we leaved the 
notion y for the function and x for the independent variable. The [y, z)-projection of 
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(^OD after a simple transformation the system can be given in the form of equations 
(m with coefficients Cly clS follows: 



1 {x + 2a)x — y 

ai =0, a2 = — , 03 = , 

Sy 6xy 

a?x — b + ay + ax"^ (2a + c)x + ba + x'^ 



xy ?/2 



(40) 



The semi invariant 1/5 has the value 

2a -\- X 

1/5 = (2ax^(18a - 2a^ - 9b + 9c) + 

9x^y^" 

x\l8a - 8a^ + 9b- 27c) + Sa'^x^ + 7ax^ + 2x^ + lOa^y^ + 9axy^) . 

It is easy to see that the invariant z/5 cannot be identically equal to zero for any 
values of the constants. Consequently, for the projection ( ^OD there does not exist 
any immersion V2{x,y) into P^(M) and the two dimensional variety V2{x,y) can be 
immeresed into P^(M) only. Also, we cannot expect to find a first integral of the 
system pD| ) in this way. We calculated also other invariants of this series, namely 
z/7 with weight 7 and z/g with weight 9. For shortness we represent here only the 
formula for z/7 (that for ug is too large): 

U7 = 273;V5 (360a^fe(-18a + 2a^ + 96 - 9c)x^ + 360a^(-36a^ + 4a^ + 

5a% - 3Qac + 2a^c + 186c - 9c^)x^ + 90a(-216a2 + 48a^ + 

18a6 + 2a% - 96^ - 36ac + 20a^c + 276c + 36c^)x^ + 

90(-108a^ + 24a^ - 18a6 - 17a^6 + 72ac + 2a^c - 96c + 27c^)x'' + 

90(-18a - 32a^ - 96 - 110^6 + 27c - 17a^c)x'^ - 90a(39a + 26 + llc)x^ + 

(-1350a - 180c)x^ - 180x^° + 10x^(2a + xf{36a^ - 4a^ - 

18a6 + 18ac + 18ax - 8a^x + 9bx - 27cx + 3a^x^ + 7ax^ + 2x^)y + 

27x{2a + x){-AOa^b - 80a^x - 3Qa%x - AOa'^cx - UAa'^x^ + lOabx^ - 

bOacx'^ - A2ax^ + 4a^x^ - 36x^ + 9cx^ + 4a^x^ + ax^)y'^ + 

14ax(2a + x)^(10a + 9x)y^ - a(200a^ + h32a^x + A8Qax^ + I3hx^)y^) . 

Also in this case we cannot reach z/7 = identically by any choice of the parameters 
a, 6, c. The existence of a term which does not containg parameters allows us to 
conclude that the other members of this series have the same property. Naturally, 
all these formulas will gratly simplify in the case a = 0. We can expect, that the 
behaviour of the system (^) in that case will show some special features. 

We take into account also the second projection of (^), i.e., that on the {x,y)- 
plane. The corresponding second order differential equation has the form (Q) and 
the coefficients ai{x,y), i = 1, 4 are given by: 

1 ax + c — a — y 

ai = 0, a2 = — , 03 = , 

3y 3y 

ax^ + cx + b (a^ + l)x + ac— 1 
04 = h a, (41) 

y 
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where a, b, c are constants of the system pl|). The calculation of the semi invariant 
z/5 results in 

— Qj ~\~ C ~\~ ddC 

z/5 = T^ (18a26-36a6c + 186c2 + (-360^6 + 18a2c + 36a6c-36ac2 + 

l%c^)x + (ISa^ + ISa^fe - na^c + bAac^)x^ + (-360^ + 54a2c)x3 + 18a=^x^ + 
(4a^ + 45a6 + 27ac - 16a^c - 366c - 21c^ + 2Aa^c^ - IQac^ + 4c^ + 
(54a^ - 16a^ - 36a6 - 36ac + 48a^c - 36c^ - 48a^c^ + lQac^)x + 
(-Qa^ + 24a^ - 72ac - 480^0 + 2Aa'^c^)x^ + (-360^ - 16a^ + lQa^c)x^ + 
4aV)?/ + (-12a^c + 24ac^ - 12c^ + (-12a^ + 48a^c - 36ac^)x + 
(24a^ - 36a2c)x^ - I2a^x^)y'^ + 

(-a^ - 8ac + 8c^ + (-80^ + lQac)x + 8aV)?/^) . (42) 
The projective curvatures Li and L2 are given by 

Li = — -(9(ax^ + cx + 6) + 2y(a — c)^ — 4ax?/(a — c) + 2a^x^y + 
3|/4 



?/^(a — 2c — 2ax)), 
c — a + ax 



(43) 



The semi invariant z/5 will be equal to zero only for a = c = with arbitrary h. In 
this case the system (^T]) has a very simple form as well as the projective curvatures 

L^ = % L, = 0. 

The calculation of the next possible series of invariants gives us using formula (|ll|) 

wi = 0. 

Consequently we must turn over to the series with the first invariant beeing ^2 



1|,|T6D. For this semi invariants the formulas hold 



. _ 156 . _ 1356^ . 24306 



y5 ' ylO ' yl5 

. _ 656106^ . 23619606^ 
^8 — ^ — ) Ho 



and the corresponding absolute invariants (|l^ are given by 

_ 3 _ 18 _ 162 _ 1944 
J4-^, J6-^, Js-Y^, W- 

Because of L2 = we cannot expect to find a first integral of the system in the form 
of (|T^. But we see, that in this case V2{x,y) can be immersed into P^'' 
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Let's mention that for a = 6 = c = both curvatures (^31) equal zero and that 
therefore our equation (41) can be reduced to y" = and trivially solved. 

Now we turn to the last case - the (x, z)-projection of the system (P0|). The 
coefficients of the corresponding equation are 

1 b — 3y — a 

ai = 0, aa = — , as = , 

3y 3ky 

a 2{b + y) x+1 cx — b — ax"^ ab 

a4 = ^ \ \ 5 • 

X x^ y y^ xy 

The semi invariants z/5 has in this case the form 

= 1^^-1^(26^ + ab'x + 9b^x^ - 6a%'x' + x\-9ab + a% - 9b' + 96c) + 

x\-9a' + 2a^ - 18a6) + ISa^x^ - 2y%'' + Safexy^) 

and the curvatures are 

— 96x^ — 9cx^ — 9ax^ + 26^?/ — Aabxy + 2a^x^y — 2by'' + 3axy^ 



Sx'^y'^ 

b — ax 



xy^ 



We also investigated the invariants z/5, 1/7, ... in the general case (a, b,c ^ 0). As an 
example we note here 

z/7 = — ^—(1806V-906^ (a + 2c)x^ + 906^(9 -7a2 + 2a6 + ac)x^ + 
27x'?/^^ 

906^ (-18a + 7a^ -9b- a% + 7a2c)x^ + 9%'{a'^ - 7a% + 9ac - 7a^c + 

9bc - 9c')x^ + 90a6(9a2 - 2a^ + 18a6 + 7a^6 - 96^ - a^c + 186c + 9c^)x^ + 

90a2(-18a6 + a% - 96^ - 9ac + 2a^c - 456c)x^ + 

90a^(9a - 2a^ + 36b + 18c)x^ - 1620a^x^° + 10(-6 + ax)^(26^ + 

ab^x + 96^x^ - Qa%''x' - 9abx^ + a^bx^ - 96^x^ + 96cx^ - 9a^x^ + 

2a^x^ - 18a6x^ + 18a'x^)y + 27x{b - axf{ab^ - 6b'x + 5abx' - 

aV - 36^x^ + 116cx^ - 20a6x^ + 2a' x%'' + 

Ub{-b + ax)^(-26 + 3ax)y^ + (-86^ + 17a6^x - 9aV)y^) . 

We can reach z/5 = only for a = b = 0. The semi invariant z/7 = only for a = b = 
in this case too as well as for all following invariants of this series. Then we have 
Li = 3cx/y^, L2 = and wi = 0. The semi invariants of the series (|15[)-(|T6]) in this 

case are 

15cx 135c2x2 2430c3x3 



^2 — , ^4: — T7\ , H 



5 ; "4 in ' " 15 ) 

y y y 

65610cV 2361960c5x^ , , 

^8 = ' ^10 = (45) 

The corresponding absolute invariants are exactly the same as in the former case 
). In the considered case V2{x,y) can be also immersed into P^'' 
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4 Conclusion 



Wc have investigated some projections of well known dynamical systems - the 
Lorenz- and the Rossler-system. The proposed new approach to the investigation 
of those systems gives us a simple possibility to distinguish the regions of the pa- 
rameters for which the systems show regular behaviour. Also, it turned out to be 
possible in some cases to find some first integrals of the considered systems. These 
observation lets us hope that the method suggested can be used as a first step for 
the investigation of further new systems. Using the proposed in this paper method 
we can distinguish different areas in the space of parameters which allow a chaotic 
behaviour. This knowledge can be used as input for different, regular or numerical 
for instance, methods. The present paper is naturally a first step towards a more 
general approach which allows to find an analytical criterium for the parameters 
of the system under what the system has a certain dynamical state, a regular or a 
stochastic one. 
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